Effect of an Electron-phonon Interaction on the One-electron 
Spectral Weight of a d-wave Superconductor 
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Here we analyze the effects of an electron-phonon interaction on the one-electron spectral weight 



A{k,iL>) of a d^ 



superconductor. We study the case of an Einstein phonon mode with various 



momentum-dependent electron-phonon couplings and compare the structure produced in A{k, uj) 
with that obtained from coupling to the magnetic 7r-resonant mode. We find that if the strength of 
the interactions are adjusted to give the same renormalization at the nodal point, the differences in 
A{k,Lj) are generally small but possibly observable near k = (tTjO). 
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I. INTRODUCTION 

The role of the electron-phonon coupling in the high 
Tc cuprates remains a puzzle. The initial finding of the 
absence of a phonon signature in the temperature depen- 
dence of the resistivity^ and the small size of the isotope 
effect in the optimally doped cuprates^ suggested that 
the electron-phonon interaction played a relatively unim- 
portant role in these strongly-correlated materials. How- 
ever, large isotope effects away from optimal doping^Sii 
significant phonon renormalization induced in the su- 
perconducting state^i^iSii and recent interpretations of 
ARPES data-*^ continue to raise questions regarding the 
nature and role of the electron-phonon interaction in the 
high Tc cuprates. 

One point of view is that the effects of the strong 
Coulomb interaction act to suppress the electron-phonon 
interaction and that while the electron-lattice interac- 
tion enters the problem, it does so on a secondary level 
coming along as it were for the ride. For example, in 
this view the large isotope effects observed in some of 
the cuprates away from optimal doping arises from the 
influence of the lattice on stripe fluctuations, acting to 
stabilize these and thus suppressing superconductivity.^" 
Similarly, the superconductivity-induced phonon renor- 
malization and the possible Englesberg-Schriefferi!- sig- 
nature in the ARPES data could be interpreted as natu- 
rally occurring in an interacting system but having little 
effect on the underlying superconducting pairing mech- 
anism. Alternatively one might interpret the isotope ef- 
fect and the phonon renormalization as supporting the 
existence of a significant electron-phonon coupling. Fur- 
thermore, ARPES measurements have been specifically 
interpreted in terms of phonon modes that could drive 
dx2-y2 pairingi^ Here, we analyze a simple model of an 
electron-phonon interaction with the goal of obtaining 
insight into what one expects to see in the ARPES data 
of a dx2-y2 superconductor with electron-phonon inter- 
actions. 

Continuing technological advances along with im- 



proved sample quality have allowed angle-resolved pho- 
toemission spectroscopy (ARPES) to probe details of the 
energy and momentum structure of the one-electron ex- 
citations in the cuprate materialsii^ Although simplified, 
the sudden approximation leads to a useful picture in 
which the ARPES intensity is equal to the square of a 
matrix element which depends upon the photon energy, 
polarization, and the sample geometry times a product 
of the single particle spectral weight 



A{k,uj) = - Im{G(fc,w)} 

TT 



(1) 



and a Fermi factor f{u!). Here G{k, uj) is the one-electron 
Green's function. Thus, the idea is that from the k and 
u! dependence of the ARPES data, one can extract infor- 
mation about the spectral weight A[k,Lu). Then, from 
this, one seeks to learn about the electron self-energy 
S(/c,w) and the structure of the effective interaction. 
In particular, the role of spin fluctuations and the tt- 
resonance on the superconducting state spectral function 
have been studiediiii^ii^ With the recent suggestions^ 
from ARPES measurements that there may be a signifi- 
cant coupling of the electrons to a phonon with an energy 
of order 40meV, one would like to understand how this 
would effect the ARPES spectrum. 

From the number of atoms in a unit cell, it is clear 
that there are a large number of phonon modes in the 
cuprates. Here we will focus on several of the modes asso- 
ciated with the motion of the O ions. We will treat these 
as Einstein phonons. Then for a Hubbard- like model in 
which the Cu sites form the Hubbard lattice, the effective 
electron-electron interaction is 



Viq,L0) 



2\giq)\^no 
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(2) 



If |g(g)P = I5P is independent of the momentum trans- 
fer, V{q,uj) does not couple to the (i2.2_j,2 -pairing chan- 
nel. This could model the coupling to the c-axis vibra- 
tion of the apical oxygen. Alternatively, if the electron- 
phonon matrix element is momentum dependent, the 
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interaction given by Eq. (0) can couple to the d^^-y^- 
pairing channel. 

The possibility that an electron-phonon interaction 
could give rise to d-wave pairing has been discussed by 
various authorsAiLiLiSii^iSSiSi In one approach, the d- 
wave pairing interaction occurs as the result of the in- 
terplay of the O half-breathing mode and the exchange 
interaction.^ Other approaches suggest that the Coulomb 
interaction can lead to a peaking of the electron-phonon 
coupling at small momentum transfers which favors 
dx2_y2 pairing^SiSS This type of momentum dependence 
also occurs directly for certain phonon modes. For exam- 
ple, for the Cu-O-Cu buckhng-hke mod9iii2J,2i^2i^ the 
square of the electron-phonon coupling constant is 



|g(<z)P = |ffP (cos2(|)-,cos2 (I)) 



(3) 



Setting q = k — k', the momentum-dependent part of this 
coupling factors into a sum of separable terms 

|5(fc-fc')P- 

+ -(cosfca; — cosfcj^)(cosfc^ — cosfc^) + ■ ■ - ^ , (4) 

including additional (cos fca; -f cos fcj, ) and (sin/caji sin ky) 
factors. The plus sign in front of the d-wave term implies 
that this type of phonon exchange provides an attractive 
channel for d-wave pairing. The key point is that if the 
electron-phonon coupling \g(k, k')\ falls off at large |fc— 
momentum transfers, then such a phonon exchange can 
mediate d-wave pairing. 

Alternatively, an in-plane O breathing-like mode has 
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('Z)P = l5p(sin^(|)+sin^(|)). (5) 

This increases at large momentum transfers giving rise 
to a repulsive interaction in the da;2_j^2-channel. Setting 
q = k — k' in Eq. ||SJ) one finds that 



\g{k-k')f = 



gp^l — - (cosfc^ — cos/cj,)(cosfc^ 



cos k', 



(6) 



and the minus sign in the second term implies that this 
phonon suppresses d-wave pairingiiLiSiiSiSSiSi 

In Section II we discuss the simplified case of a cylin- 
drical Fermi surface and a separable phonon mediated 
interaction. This provides insight into the differences be- 
tween the s-wave and d-wave cases and establishes the 
structure of the singularities in the self-energy that are 
reflected in A{k,uj) for an Einstein mode. While in the 
actual materials, the singularities are broadened by the 
dispersion of the phonon mode, quasiparticle lifetime ef- 
fects due to other interactions and impurities, as well as 
finite temperature effects, these results provide a useful 
framework for understanding the structure that appears 
in A(k, Lu). 

In Section III, we include the effects oi a t — t' band 
structure and the momentum dependence of the cou- 
pling. We consider the three different electron-phonon 



coupling constants discussed above and compare these 
with the response to the 7r-resonance spin fluctuation 
mode. The analysis of the 7r-resonance mode has been ex- 
tensively discussed in Refs 14"'15. Various estimates for 
the strength of the 7r-resonance-electron coupling have 
been madei2i24 Here, we chose this coupling so that the 
renormalization of the nodal Fermi velocity is compara- 
ble with that obtained for the phonon coupling. Then we 
compare and discuss the spectral weights for the various 
modes. Section IV contains a summary of the results and 
our conclusions. 



II. A CYLINDRICAL FERMI SURFACE AND 
AN EINSTEIN PHONON 

In this section we consider the case of a cylindrical 
Fermi surface and an interaction arising from the ex- 
change of an Einstein phonon of frequency f2o 



2\g{e,e')\^no 



(7) 



Here, 9 and 9' denote different k vectors on the cylindri- 
cal Fermi surface. With Eq. Q in mind, we will take 
\g{9,9')\'^ to have the separable form 

\gi9,9')\^ = \g,f + 1 5^ p cos 20 cos 20'. (8) 

The one-electron Green's function can be written as 

uj + tk 



G{k,uj) 



{Ziuj)Lor-el-^^9,Lu)-' 



(9) 



with efe — fc^/2m — ^, Z{uj) the renormalization parame- 
ter and (p {9, lo) = cj) [uo) cos(20) the gap parameter. The 
Eliashberg equations for Z{ljj) and {lj) are 

(l-ZH)c. = A,/„°°da-'/f 



Re 



[(Z(w') w')2-02(j^/)(.os2 29] 2 

( I u 

\ a;'+cj+£lo— i (5 lj' —u)-\-ilQ—i 5 J J 



(lOa) 



0H = A,/„°°da;'/f 



Re 



0(m')cos^ 20 



( ^ 
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uj' —LU-\-il.(j —i 5 



)} 



(lOb) 



with K = 2|5,|2iV(0)/f7o and = 2\g^\'^N{0)IQ.a. Here 
A^(0) is the one-electron density of states at the Fermi 
surface. 

In order to determine the effect of the phonon on Z{io) 
and we will adapt an approximation used in the 

early studies of the role of phonons on the superconduct- 
ing I(V) characteristiC(2^ From the form of eqs HlUa|) and 
(|10b|l . one sees that there will be structure in Z{uj) and 



(j) {u) when w ~ ± (JIq + A(0)). In this case, uj' will be of 
order the gap A{9) at the gap edge 



5.0 



(f){0,uj = A{e)) 

Z{6,u)^ A{0))' 



(11) 



Therefore, if the low-energy response in the supercon- 
ducting state is well described in terms of BCS d-wave 
quasiparticles, one can replace Z{uj') and (j){9,LLj')/Z{uj') 
inside the integrals by Z{0) and A{9) = Ao cos 20. Then, 
taking the imaginary parts of eqs H10a|l and IjlObjl we have 
for cj > 



ujZ2{uj) = AXz 
<f)2{i^) — 4Aa 



dO- 



- (12a) 



e. [{uj-VLoY - Alcos'^29] 

Ao cos^ 20 , „, , 

d9 (12b) 

Sc [{uj -n^f ~ Alcofi'^29\^ 



Here 9c is such that A{9c) — uj — and (j)2{uj) and 
^2(0;) are even functions of cu for a time-ordered zero 
temperature Green's function. 
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FIG. 1: Results for the real and imaginary phonon-induced 
contribution to Z{uj) and (j}{io) for a d-wave (solid) and an 
s-wave (dashed) superconductor. Here we have taken a cylin- 
drical Fermi surface and a separable interaction. Z{u>) and 
<j!>(cj) are normalized with respect to the appropriate coupling 
constant, X,/, and Xz for the d-wave gap and Xz for the s-wave 
case. For all of the circular Fermi surface plots, energy is 
measured in units of Ao and fio = 1-5 Ao. With this normal- 
ization 0^^2(0;) goes to 7r/2 as a; ^ cxd. 



Results for ^^2(0;) and 4>2{^) are shown in the top 
panel of Fig. 1 for both a d2,2_y2-wave and an s-wave 
gap with r^o = l-5Ao. For an s-wave gap, cos 20 is set 




FIG. 2: An intensity plot of the spectral weight A{k,uj) for 
the case of an s-wave superconductor coupled to a phonon 
with frequency f2o- As indicated on the color scale, the figure 
is 'overexposed', i.e. features exceeding 10% of the maximum 
intensity appear white, in order to show the weaker features. 



to 1 and = in eqs H12a() and (|12b|) . For the s-wave 
case, the imaginary parts of Z{ijj) and 4>{ijj) onset when 
UJ exceeds ± (fio + Ao) and exhibit a square root singu- 
larity. For a d2;2_j^2-gap, these functions onset linearly at 
= ± i7o because of the gap nodes and there is a log 
singularity at ± (Slo + Aq). The real parts of Z{uj) and 
(j) {uj) are obtained from the usual dispersion relations, 
and results for Zi{uj) and (/>i(w) are shown in the lower 
panel of Fig. 1. For the s-wave case, 0i and Zi exhibit 
square root singularities as uj approaches ± {VLq + Ao). 
This is just the expected Kramers-Kronig transform of 
the square root singularity in 02 and Z2. Similarly, the 
results for (j)i and Zi for the d^2_y2 case exhibit step 
discontinuities at w = ± {Hq + Aq) arising from the 
log singularities in 02 and Z2. Naturally in real mate- 
rials, phonon dispersion, impurity scattering, and finite 
temperature effects broaden these features. Neverthe- 
less, they provide a simple framework for analyzing the 
ARPES data. 

An intensity plot of A{k,uj) for the case of an s-wave 
gap is shown in Fig. 2. Here, A{k,uj) is obtained from 
the imaginary part of G{k,uj), using the s-wave results 
for Z{uj) and (j>{uj) shown in Fig. 1 with A = 0.5. The 
real part of the gap function is supplemented by an ad- 
ditional contribution from an underlying pairing interac- 
tion so that the magnitude of the gap at the gap edge 
is equal to Aq. Results for both the ARPES accessi- 
ble region uj < and the inverse photoemission region 
a; > are shown. The shift of spectral weight due to 
the quasiparticle coherence factors ^(1 -I- ^) is clearly 
seen as is the Englesberg-Schrieffer signature showing 
the asymptotic approach of a peak in the spectral func- 
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FIG. 3: Energy distribution curves show 
u) for various values of for an s-wave sr 
Az = 0.5. Here one sees that as 6^ exceed 
peak is left behind whose intensity weaker 



tion to ± (flo + Ao). Because of the 
gularity in Z and 0, the asymptotic 
peak to ± (ilo + Aq) varies as (Ailo/efc 
the Fermi velocity is renormalized by 
that the dispersion of the peak for w : 
A/(efc/2'i(Ao))2 + while for w large 
a broadened quasiparticle peak disperf 
distribution curves (EDC) showing 
various values of are shown in Fig. : 
is the type of EDC that one would e> 
traditional s-wave electron-phonon sup^ 
a single dominant Einstein mode.^^ More generally, one 
would have multiple phonon modes and their dispersion 
along with possible finite temperature effects would lead 
to a richer response. 

Intensity plots of A(k,Lo) for the case of a dx-2_yi 
gap are shown in Fig. 4. Just as for the s-wave case, 
01 (0, lS) is supplemented so that the gap at the gap edge 
is Aocos26'. Fig. 4(a) shows Aik^uS) for a cut along the 
antinodal direction in fc-space (Q = 0), while Fig. 4(b) 
shows the results for a cut along the nodal direction 
{6 = 7r/4). The antinodal cut resembles the s-wave case 
in the transfer of spectral weight as passes through the 
Fermi energy and the renormalization of the quasiparti- 
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FIG. 4: An intensity plot of A{k, uj) for a d-wave gap. Results 
for a. 6 = cut are shown in (a) and in (b) for 9 — it /A, a 
nodal cut. 



cle dispersion. However, the Englesberg-Schrieffer signa- 
ture no longer asymptotically approaches ± (fio + Aq), 
but rather appears to be broadened and cut off. In the 
s-wave case, the broadening due to the electron-phonon 
interaction did not set in until exceeded f2o-f- Aq lead- 
ing to the long sweep of the peak which occurs for |cl;| just 
below {flo + Aq). However, the nodal regions associated 
with a dx^-y'i gap lead to a finite broadening when |a;| 
exceeds r^o. The onset of this broadening is seen by the 
faint horizontal line in Fig. 4 where the intensity changes 
from black to blue at larger values of efc. As we will dis- 
cuss, termination of this peak is a reflection of the fact 
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FIG. 5: Energy distribution curves (EDC) showing A{k, u) 
versus uj for different values of ep with Xz — X4, = 0.5. Results 
for the 9 — cut are shown as the solid curves and for 9 — tv/A 
as the dashed curves. 



that for a dx^-y^-gs-p, Zi and (f>i have step discontinuities 
at ± (f2o + ^0) rather than the square root singularities 
associated with an s-wave gap. 

The nodal cut, shown in Fig. 4(b), appears on first 
glance to be similar to what one would expect for the 
normal state. That is, a renormalized ek/Zi{k,0) disper- 
sion for w <C rio with the dispersion returning to its band 
value efc for oj ^ ^Iq. However, the cut-off Englesberg- 
Schrieffer signature still occurs for \uj\ ~ flo + Aq. Thus, 
the full antinodal gap Aq enters as the characteristic kink 
energy for all momentum slices. This simply reflects the 
w| = r^o + Aq singularities in Z and 4> shown in Fig. 1. 
Again, the broadening of the Englesberg-Schrieffer peak 
when \uj\ exceeds Hq is clearly seen in Fig. 4(b). In Fig. 5, 
various EDC slices of A{k,uj) are shown for the dx2_y2 
case. Comparing these with the s-wave case, one sees the 
broadening and truncation of the Englesberg-Schrieffer 
lower peak. 

The difference in the structure of the Englesberg- 
Schrieffer signature between the s- and the da,2_j,2 -cases 
can be understood from the plots of 

Ek = -^{Z,{ij)u;)^ -cf,l{u;) (13) 
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FIG. 6: Plots of £fc = —\/ {Zi{uj)uj)'^ — ^) versus lu along 

the negative uj axis for = 0.5 show the structure of the 
Englesberg-Schrieffer signature. The upper panel is for an s- 
wave gap and the lower panel is for a d-wave gap with 9 = 
shown as the solid curve and ^ = 7r/4 as the dashed curve. 



shown in Fig. 6. One can see that as one probes Ck states 
which are further below the Fermi energy, two solutions 
of Eq. p3|l develops. For the s-wave case shown in the 
upper panel of Fig. 6, an undamped lower energy branch 
asymptotically approaches lo = — (Slo + ^o); ^ind a second 
quasiparticle branch at w ~ — evolves which is damped 
by the imaginary parts of Z and (j). As we have seen, 
these branches are reflected in the structure of A{k,uj) 
and the lower energy branch represents the characteristic 
Englesberg-Schrieffer signature for an s-wave supercon- 
ductor. Similar plots for the d2.2_j,2-case with 9 — and 
6 = tt/A are shown in the lower panel of Fig. 6. Here, 
unlike the s-wave case, the low energy branch is termi- 
nated, reflecting the fact that the singularities in Zi and 
4>i for the d-wave case are simply step discontinuities at 
± (fio + Aq). The onset of damping processes for the 
dj.2_y2 case when u < — Oq give rise to the discontinuity 
in slope seen at uj = —Qq. 



III. BAND STRUCTURE AND THE EFFECT 
OF A MOMENTUM-DEPENDENT COUPLING 

We turn next to the effects of the band structure and to 
the momentum dependence of the electron-phonon cou- 
pling. For the band structure, consider a square lattice 
with a near-neighbor hopping t and a next-near-neighbor 
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hopping t' . In this case 

€k — — 2i(cos kx + cos ky) — 4i' cos k^ cos ky — fi. (14) 

For t' /t — —0.3 and fi/t — —1, one has the typical Fermi 
surface shown in Fig. 7 and the single spin electron den- 
sity of states shown in the inset. We take the gap to be 

Afc = Ao(cosfc2. — cosA:j^)/2. (15) 

In this case, the one-electron Green's function can be 
written in the form 

^ Z{k,uj)u; + {ek+X{k,uj)) 

^ '""^ iZik,u;)Lor-{ek+X{k,u;))^-<j)^k,iuy 

(16) 

Adopting the same approximation as before, the phonon- 
induced contributions to the imaginary parts of the 
renormalization, energy shift, and gap parameters are 
given by 

coZ^ik^Lu) = -^^|g(fc-fc')P X 

k' 

{s{Ek' + no-u;)-S{E„, +no+u;)y (17a) 




k/7t 

FIG. 7: The Fermi surface for t' /t = -0.3 and = -1. The 
inset shows the single spin electron density of states. We will 
discuss the spectral weight for the cuts marked A, B, and C. 



X,ik,u;)^-^ Y^\g{k^k')\' 

((5 {Ek' +na-uj) + S {Ek' +no+ 1^)) , (17b) 



Mk,io) = - Y,\gik-k')\' 




X 



(,5 [Ek' +Qo-i^) + S (Ek' +no + w)) • (17c) 

Here, as before, we assume that an underlying pairing 
interaction, most likely spin-fluctuations, gives rise to 
a zero temperature d.j.2_y2 superconducting state. At 
low energies this state is characterized by a renormal- 
ized band structure Eq. H14|l and chemical potential, a 
renormalized coupling constant g{q), and a gap given by 
Eq. IjlSI) . These parameters have been used in the Eliash- 
berg equations to describe the state which enters when an 
excitation at energy uj > fig decays to a lower energy Ek' 
state (or when uj < — fio decays to —Ek')- The real parts 
of Z{k,u}), (j){k,uj), and X{k,uj) are again found from the 
Kramer-Kronig dispersion relation. The spectral weight 
A{k, uj) is then obtained from Eq. ^ with the chemical 
potential shift removed from Xi{k, lu) and a contribution 
added to the real part of the gap so that the real part of 
the gap at the gap edge remains equal to A^, Eq. (fT3|l . 
Note that the contributions of the underlying pairing in- 
teraction to Z and X, as well as the higher energy part 
of 4>, have not been included. Thus, there are additional 
renormalization and damping effects which do not ap- 
pear. We basically are seeking to understand the leading 



contribution of the electron-phonon interaction which is 
superimposed on top of the other many-body interac- 
tions. This approach rests on the idea that in the super- 
conducting state the low-lying electronic states are well 
described by BCS dx^-y2 excitations^^ with renormalized 
band parameters t, t', and /i, a dx2-y2-wave gap Ak, and 
renormalized electron-phonon coupling constants. Note, 
that here we are not taking into account the possible 
change in g-dependence of the electron-phonon couplings 
produced for example by the Hubbard U jiLS^ 

We begin by looking at the self-energy terms for the 
case of the buckling mode with |g((z)|^ given by Eq. Q 
and |(7p = 0.5 in units of t^^. Results for (l){uj,k), 
Z(uj, fc), and X(uj, k) are shown in Fig. 8 for k at point A 
shown in Fig. 7. The imaginary parts of Z and </> exhibit 
the expected log singularity at Ao + r2o that we previ- 
ously saw for the case of a circular Fermi surface. In 
addition, there is a second log singularity at -E(0, tt) -|- fio 
with E{Q, tt) = v/e^(0,7r) A2 which comes from the 
Van Hove singularity^^ at A; = (0,7r). These log singular- 
ities in Z2 and 02 manifest themselves via the Kramers- 
Kronig dispersion relation as step-down discontinuities in 
Zi and (/ii, as seen in Fig. 8. The energy shift parame- 
ter X has only the Van Hove singularity. Naturally, the 
dispersion of the phonon mode as well as finite temper- 
ature and lifetime effects will broaden these features in 
the actual system. The energy distribution of the spec- 
tral weight A{k,ijj) for the buckling mode at momentum 
kA is plotted in Fig. 9. It shows the quasiparticle peak 
at the gap edge A^^ as well as structure associated with 
the buckling phonon at fio + Ag and fio + £^(0, tt). 

As noted in the introduction, one would like to de- 



7 




FIG. 8: The self-energy parameters (j), Z, and X versus uj 1 
the case of the buckhng phonon couphng, Eq. Q , with k 
\^A corresponding to the point A of Fig. 7. Here = 0.5 



termine whether the structure observed in the ARPI 
data is due to phonons or the 7r-resonance spin fluctu 
tion mode. Eschrig and Normaniiii^ have analyzed t 
effect of the 7r-resonance using a detailed tight binding 
of the band energy Cfe and a coupling to the 7r-resona 
mode of frequency f2o given by 



^^4^2[cos2(g^/2)+cos2(g^/2)]- ' 

Here, we will use the t — t' band structure of Eq. Ijl I 
with t' /t — —0.3 and /i = —1, set wq = 1, ^ = 2, aj 
set g^-p = 5 which corresponds to having = 
in an effective Hubbard RPA interaction. In additio 
with this choice for ggp we will find that Zi(fcF,0) 
the nodal point C is comparable with Zi{kF,Q) for t 
phonons. This makes it convenient for addressing t 
question of whether there are significant spectral d 
ferences due simply to the structure of the momentui 
dependent couplings that would allow one to determi: 
the nature of the mode from the ARPES data.^^ No 
that for the spin-fluctuation interaction with giv 
by Eq. (|18|l . there is a minus sign on the righthand si 
of Eq. (|17c|l for the gap parameter. For the three types 
phonon couplings we take \g\'^ = 0.5 in units of i^^. Tl 
gives Zi(fcF,0) ~ 1.3 corresponding to an effectiv o^"- 
A ~ 0.3. For the 7r-resonance mode coupling, setting 
alp = 5 gives Zi{kF,0) c^i 1.3. 

Intensity plots of A{k, ut) for the constant Holstein cou- 
pling, the buckling mode coupling Eq. the breathing 




FIG. 10: Intensity plots of A(k,u)) for the momentum 
cut A for four different couplings corresponding to the Hol- 
stein mode with \g{q)\'^ constant, the buckling mode Eq. lO, 
breathing mode Eq. JSJ, and vr-resonance magnetic mode cou- 
pling Eq. (|THJ. Here, \g\^ = 0.5 and Qsf = 5. The cut-off in- 
dicated on the color scale refers to the actual spectral weight 
intensity (as opposed to the relative scale used in the previous 
intensity plots) so that one can directly compare the effects 
of the different couplings. 




0.0 kx/ic 33 0.0 [^^/x 0.33 



A[k,a] 3J) 



FIG. 11: Intensity plots of A{k,u!) for the momentum cut 
for the four different couphngs. 



mode coupling Eq. and the 7r-resonance mode co 
pling Eq. (0, are shown in Fig. 10 for the momentum 
cut A. Similar intensity plots for the momentum cuts 
B and C are shown in Figs 11 and 12. In Fig. 10, one 
sees a high intensity quasiparticle peak and weaker struc- 
tures onsetting at uj = — (JIq + ^o) and — (ilo + £'(0,7r) 
due to the couphng to the phonon or magnetic resonance 
modes. For the B momentum cut shown in Fig. 11, one 
can now move deep enough inside the Fermi sea that 
the Englesberg-Schrieffer lower energy peak (the upper 
bright curve in the figures) is broadened when lo becomes 
less than —Qq and terminated at a finite value of kx as to 
approaches — (fio + Aq). At still higher energies (w more 
negative), a damped quasiparticle branch is seen. The 
nodal C cut is shown in Fig. 12. Here, one clearly sees 
the Englesberg-Schrieffer signature with a quasiparticle 
peak which varies as ek/Zi{kp, 0) near the Fermi surface, 
then disperses and bends as oj approaches — (f^o + Aq). 
This peak is then terminated as a broadened high energy 
quasiparticle branch appears at more negative values of 

LJ. 

The difference of A{k,uj) for the various modes is in 
fact subtle since all four have an Einstein spectrum with 
ilo = 0.3t, a dg.'^_y2 gap with Aq — 0.2t, and a band 
structure with t' /t = —0.3 and fi — —1. Thus, the 
characteristic energies Aq, flo + Aq, and flo + E{0,tt) 
are the same. In addition as discussed, we have chosen 
the couphng constants so that |5(g)P averaged over the 
Brilloin zone is the same for all four cases. Thus, the 
basic difference is the momentum structure of the dif- 



Holstein phonon Buckling phonon 
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FIG. 12: Intensity plots of A{k,uj) for the nodal momentum 
cut C for the four different couplings. 




qy7l=q/7l 

FIG. 13: The momentum dependence of the |g(?)P couphng 
versus q = qx = Qy for the four different modes with \g'^ = 0.5 
for the phonon modes and ggp = 5 for the vr-mode. Note that 
this is a slice of a two-dimensional {qx,qy) surface and that 
the volume enclosed by these surfaces is {2'K)^\g\^ for each of 
these couplings. 
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FIG. 14: Energy distribution curves (EDC) showing A{k, u) 
versus tj at fc = (0, vr) for the four different modes. 



ferent couplings shown in Fig. 13 for = Qy Here, 
we see that the spin-fluctuation resonant mode is clearly 
most strongly peaked at large momentum, followed by 
the breathing mode phonon, the uniform Holstein cou- 
pling, and lastly the buckling mode phonon which has 
|g((7r, 7r)p = 0. One consequence of the strong peak in 
the magnetic reasonance-mode coupling is seen in Fig. 10 
for the A cut. Here, the increase of the intensity of the 
spectral weight A{k,uj) which occurs when uj decreases 
below — (rio + E{tt, 0)) is greatest for the spin- fluctuation 
TT-resonance and smallest for the buckling mode. 

In Fig. 14 we show the energy distribution curves for 
the four modes for momentum k = (0,7r). for 
all of the modes shows a strong peak at Aq. For the 
TT-mode, this is followed by a dip and then a secondary 
peak which develops as w decreases below the Van Hove 
threshold at — (fio + ^'(0,7r)). It is this peak-dip-hump 
structure, for the case in which the effects of the bilaycr 
splitting can be eliminated, that has been identified as 
a 'fingerprint' of the 7r-resonanceii^4i^*^2i2^ Here, we see 
that indeed this structure is most pronounced for the 
TT-mode and smallest for the buckling mode. However, 
this is a quantitative effect rather than a qualitative one 
and if the phonon coupling increases at large momentum 
transfers, such as in the case of the breathing mode, this 
feature returns although not as strongly as for the tt- 
mode. 



IV. CONCLUSIONS 

The Englesberg-Schrieffer-like structure in the ARPES 
data of BISCO is consistent with the existence of an 
Einstein-like mode with ^Iq ^ 40meV coupled to the elec- 
trons as suggested by various authors However, it 
seems that it will be difficult to determine the origin of 
the mode based solely upon the {qx,qy) momentum de- 
pendence of its coupling. One might have thought that 
the g-dependence of the coupling or in the case of the tt 
mode, the g-dependence of the resonance that has been 
parameterized as a g-dependent coupling, would give rise 
to clearly identifiable structure in A{k^uj). However, all 
of the modes show quite similar characteristic features at 
energies r2o + Aq and VIq+E{Q, tt), which appear through- 
out the zone. 

It would appear that the best place to look for a feature 
that could distinguish between, for example, the buck- 
ling phonon mode and the 7r-resonant mode is near the 
k = (0, tt) point. Here, the strong coupling of the tt mode 
to the electrons for q near (tt, tt) leads to a secondary 
peak onsetting at an energy u) — — (£'(0,7r) + Qo). For 
the buckling phonon mode, the coupling at q — (tt, tt) 
vanishes and there is only a relatively weak response in 
this same frequency range. However, as we have seen, 
there is a secondary peak for the breathing mode which 
has nearly the same strength as that for the 7r-mode. 
Thus, the observed peak-dip-hump structure could also 
be consistent with a coupling to the oxygen-breathing 
mode. Recently, it has been suggested that the qz de- 
pendence for a bilayer system may identify the mode as 
having qz = tt, which would provide support for the tt 
resonancCi^^ However, further work on the odd and even 
bilayer phonon coupling is needed for comparison. 

While the coupling to the 7r-resonance mode along 
with a higher energy continuum spin-fluctuation spec- 
trum provides an attractive unified framework, our re- 
sults leave open the possibility that an oxygen phonon 
mode could also play a role. As we have seen, even with 
a relatively modest coupling constant A ~ 0.3, one would 
expect to see evidence of some oxygen phonon modes. 
If they are not seen, then this suggests that the strong 
Coulomb many-body effects act to suppress the electron- 
phonon coupling. Alternatively, if it can be shown that 
the TT-mode is not viable, oxygen phonon modes could 
provide a source for the resonant mode features. The con- 
tinuum spin fluctuations would, of course, also contribute 
in this mixed scenerio. Here we should note that even if 
the mode were idcntiflcd as the buckling mode, wc find 
that its contribution to the magnitude of the (1^2 _y2 gap 
is negligible because the increase in Zi more than offsets 
the increase in (pi (in Eq. ^ the c?-wave coupling term 
is only 1/4 of the uniform coupling). This is in agree- 
ment with the results of Eliashberg-like Tc calculations^ 
which find that, while the buckling phonons can provide 
an attraction in the (i2;2_j^2 -channel, its contribution to 
Z leads to an overall reduction in Tc- 

To conclude, from the results that we have found, it 
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seems likely that the identification of the excitation re- 
sponsible for the structure in the ARPES data will be 
decided on grounds other than the momentum depen- 
dence of the effective coupling. One aspect that remains 
under discussion is the strength of the various couplings. 
For the O phonon modes, LDA calculations^L^i find in- 
termediate coupling strengths with A ~ 0.3 to 0.5. From 
our calculations it would appear that at this strength, 
one should in fact see structure in A{k,uj). If this is not 
seen, it raises the question of why is the electron-phonon 
coupling weakened in strongly-correlated materials 
The coupling to the 7r-resonance mode would appear to 
raise the opposite problem. That is, if the 7r-resonance 
mode is responsible for the Englesberg-Schrieffer-like sig- 
nature in the ARPES spectrum, how can it be coupled 



so strongly ?iiSi2^ 
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